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' Abstract. We present observations of the famous OH/IR star OH26.5-I-0.6 obtained using the Mid-Infrared 

OO . Interferometric Instrument MIDI at the European Southern Observatory (ESO) Very Large Telescope 

Interferometer VLTI. The emission of the dusty envelope, spectrally dispersed at a resolution of 30 from 8 /im 
to 13.5 fim, appears resolved by a single dish UT telescope. In particular the angular diameter increases strongly 
' ■ within the silicate absorption band. Moreover an acquisition image taken at 8.7 fim exhibits, after deconvolution, 

a strong asymmetry. The axis ratio is 0.75±0.07 with the FWHM of the major and minor axis which are 286 mas 
and 214 mas respectively. The measured PA angle, 95° ±6° is reminiscent of the asymmetry in the OH maser 

■ emission detected at 1612MHz by Bowers & Johnston (1990) for this star. In interferometric mode the UT1-UT3 
102m baseline was employed to detect the presence of the star. No fringes have been found with a detection 
threshold estimated to be of the order of 1% of the total flux of the source, i.e. 5-8 Jy. These observations were 

5^ ' carried out during the phase of maximum luminosity of the star, when the dust shell is more diluted and therefore 

^ ' the chance to detect the central source maximized. We modeled the dusty environment based on the work of 

. Justannont et al. (1996). In particular, the failure to detect fringes provides strong constraints on the opacities in 

the inner regions of the dust shell or in the close vicinity of the star. 
•i-H . 

■ Key words, radiative transfer - Techniques: interferometric - stars: AGB and post-AGB - stars: circumstellar 
' matter - stars: individual:OH26.5+0.6 



1. Introduction 

The short transition phase between the end of the 
Asymptotic Giant Branch (AGB) phase and the formation 
of a White Dwarf (WD) surrounded by a Planetary Nebula 
(PN) is still poorly understood. The drastic changes ob- 
served in the circumstellar environment of AGB and post- 
AGB stars are particularly puzzling. During the late AGB 
or early post-AGB evolutionary stages, the geometry of 
the circumstellar material of the vast majority of stars, 
changes from more or less spherical to axially symmetric, 
as shown by the large number of axisymmetric proto-PNe 
(e.g. Sahai 2000). As a result, most PNe exhibit axisym- 
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* Based on observations made with the Very Large Telescope 
Interferometer at Paranal Observatory 



metric structures, ranging from elliptical to bipolar, often 
with an equatorial waist and (sometimes multiple) jets 
(Corradi & Schwarz 1995). It is thought that pure hydro- 
dynamical coUimation provided by dense equatorial disks 
or tori (Icke et al. 1989), and/or magneto-hydrodynamical 
coUimation (Chevalier & Luo 1994) can explain the de- 
velopment of the extreme bipolar geometries observed. 
Whether these equatorial structures can arise in a sin- 
gle star scenario is still strongly debated (Bujarrabal et 
al. 2000). 

In recent years the advent of infrared spectroscopy has 
improved our understanding of the AGB and post-AGB 
evolutionary phase. The IRAS and the ISO infrared tele- 
scopes have detected large amounts of dust grains around 
these stars. The majority of observed circumstellar en- 
vironments show either an oxygen-rich chemistry or a 
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carbon- rich one. AGB/post-AGB stars dominated by O- 
rich dust chemistry are those where the third dredge-up 
never raised the C/0 ratio above 1. During the AGB 
phase, mass loss drives the evolution of the star. In par- 
ticular, mass loss increases dramatically (by a factor 10 
at least) towards the tip of the AGB in what is called a 
superwind (Iben & Renzini 1983) which ejects most of the 
remaining envelope of the star. Some 0-rich AGB stars 
exhibit OH maser emission and are called OH/IR stars 
(Wilson & Barrett, 1972). Their mass-loss rates are so high 
that the dust shell completely obscures the central star, 
and the object is observable only at infrared wavelengths 
and through molecular line emission at radio wavelengths. 
The nature and geometry of the superwind is still to be 
settled. The geometry of the maser emission is usually well 
constrained due to the combination of the spatial resolu- 
tion provided by interferometric techniques and the large 
extension of the maser (usually a few arcsec). The obser- 
vations of the youngest (i.e. more optically obscured) pre- 
planetary nebulae (PPN) where the superwind has just 
ceased suggest that asymmetries are already present. An 
extensive discussion on the appearance of bipolar outflows 
in OH/IR stars can be found in Zijlstra et al. 2001. 

OH 26.5+0.6 (RAFGL 2205, IRAS 18348-0526) is an 
extreme OH/IR star showing a large dust column density 
and hence a very high dust mass loss rate. It is one of 
the brightest OH maser emitters (Baud 1981; te Lintel 
Hekkert et al. 1989; Bowers & Johnston 1990) with a wind 
terminal velocity of 15 km s~^. Bowers & Johnston (1990) 
mapped the OH maser around the star and found a shell 
radius of about 2-3 arcsec. OH 26.5-1-0.6 exhibits a low CO 
J = 1 — and J = 2 — 1 emission (Heske et al. 1990) while 
at 10 fim, the silicate absorption indicates a large dust 
column density and hence a very high dust mass loss rate. 
The 10 fim complex is dominated by amorphous silicate 
absorption which has been studied by numerous authors. 
The ISO spectrum of OH26. 5-1-0.6 has been discussed by 
Sylvester et al. (1999) and Molster et al. (2002) studied 
particularly the signature of the crystalline silicates. 

Justtanont et al. (1994, 1996, hereafter JU96) sug- 
gested that this star has recently undergone the superwind 
phase and shows evidence of two mass-loss regimes: a su- 
perwind phase in which the mass-loss rate is 10~*MQ/yr 
which started recently (t <150 yr), and an earlier AGB 
phase with a mass- loss rate of about 10~^Mo/yr. The in- 
tegrated mass lost during the superwind phase has been 
estimated to be 0.1 Mq. 

Fong, Justtanont, Meixner, Campbell (2002) reported 
millinietric CO observations which did not show any sig- 
nificant deviation from spherical symmetry for the enve- 
lope of OH26.5-I-0.6. Nevertheless, it must be pointed out 
that the source is mainly unresolved at this wavelength. 
In contrast, it is one of the brightest and most asymmetric 
OH maser sources known among AGB stars, with a prefer- 
ential axis of symmetry oriented approximately east-west 
(Baud 1981; Bowers & Johnston 1990). 

The duration of the superwind phase depends on the 
mass that the star has to lose before the envelope is small 



enough to sustain the mechanism of stellar (photospheric) 
pulsations. For a star like OH 26.5-1-0.6, JU96 state that 
the superwind began very recently i.e. less than 150 yr ago. 
Radio emission in molecular lines is expected to change 
less rapidly than the infrared emission at the advent of the 
superwind phase. It is therefore of particular interest to 
study the mid-IR spatial geometry of OH/IR stars in order 
to determine the onset of asymmetries in the environment 
of evolved stars. 

Unfortunately, the dusty environment of OH/IR stars 
is difficult to resolve by single dish telescopes in the IR. 
Even more complicating is the time variability of the 
OH/IR envelopes which modulate their size and luminos- 
ity. OH26. 5-1-0.6 is a long period pulsating star whose pe- 
riod has been refined recently by Suh and Kim (2002) 
to P=1559±7 days. Taken into account the large vari- 
ations of the mid-IR flux from this star throughout its 
pulsation cycle, the published data on its spatial extent 
have to be systematically placed in their temporal con- 
text. Infrared speckle interferometry has been performed 
by Fix & Cobb (1988) close to the maximum. They pro- 
vide an extension for the circumstellar dust shell at 9.7 
/im (within the strong silicate absorption) at maximum of 
0.5" ±0.02", while outside this feature (at 8 /im) the shell 
remained unresolved by their experiment (at most 0.2"). 
They have also resolved the environment using the broad 
N band filter near phase 0.6 with a detected FWHM of 
about 0.3" (Cobb & Fix (1987). Some asymmetries have 
been reported by Mariotti et al. (1982), Dyck et al. (1984), 
Cobb & Fix (1987) and Fix & Cobb (1988), Starck et al. 
(1994). However, the reported asymmetries are within the 
estimated error bars of the measurements and all together 
the results are somewhat inconclusive, and sometimes con- 
tradictory. 

The Mid-Infrared Interferometric Instrument MIDI at- 
tached to the Very Large Telescope Interferometer (VLTI) 
is able to provide a spatial resolution in the mid-infrared 
ranging from the one provided by single-dish 8m telescope 
(about 300 mas) to the one provided by interferometric 
technique (about 5-10 mas). MIDI can also disperse the 
light with a spectral resolution of 30 through the entire 
N band which makes it a unique instrument particularly 
adapted for the study of dusty environments. We used the 
102 m baseline between the telescopes Antu (UTl) and 
Mehpal (UT3) to observe OH26. 5-1-0.6 for the first time. 

In Scction[51we describe the observations and the data 
reduction procedures, divided into three parts: (i) the sin- 
gle dish acquisition images fSect. lTTIl . (ii) the spatial and 
spectral information of the spectra fSect. lT^ . and (iii) the 
interferometric signal (Sect.[^J. In Sect|21we model the 
observations using a spherically symmetric dust model. 
Finally, in Section we discuss the results of our model 
fitting. 

2. Observations and data reduction 

OH26.5+0.6 was observed with MIDI (Leinert, Graser 
et al. 2003a and 2003b), the mid-infrared recombiner of 
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Table 1. Journal of observations: acquisition images 
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OH26.5+0.6 


start 


06:56:24 
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lOOs 


OH26.5+0.6 


star2 
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5000 


50s 


OH26.5+0.6 


starS 
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150s 


OH26.5+0.6 


star4 
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150s 


HD 177716 


PSF7 
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20s 


HD 177716 


PSF8 


08:04:23 


2000 


20s 



the VLTI. The VLTI/MIDI interferometer operates as a 
classical Michelson stellar interferometer to combine the 
mid-IR hght (N band, 7.5 - 14 ^m) from two VLT Unit 
Telescopes (UTs). The observations presented here vifcre 
conducted in the night of the 14th of June 2003, the UTl 
and the UTS telescopes were used, separated by 102 m 
with the baseline oriented 40° (East of North). 

The observing sequence, described extensively in 
Przygodda et al. (2003), is summarized hereafter. The im- 
ages have been recorded using the MIDI star acquisition 
modes called Default_Chop and Acquisition_chop with the 
8.7 /im filter. The Acquisition_Chop mode is the first tem- 
plate used after the pointing to test if the target is within 
the MIDI Field Of View (FOV) (diameter of about 3") and 
to perform a fine pointing. The chopping mode (f=2IIz, 
angle -90 degree) is used to visualize the star, which is 
not perfectly centered in the first image, and centered in 
a second step. The number of frames recorded per image 
is generally about 2000 and the exposure time is by de- 
fault 4 ms in order to avoid background saturation. If the 
result of the template is not satisfactory, the procedure 
is started again. It must be pointed out that no nodding 
sequences are performed, the sky being removed by chop- 
ping only. For some stars for which the coordinates are 
not well-defined, it might be difficult to get the star di- 
rectly in the MIDI FOV at the first attempt. This was 
the case for OH26.5+0.6. Therefore, the Default_Chop is 
used instead. In this mode, images are recorded only for 
visualization and the pointing is done by 'hand' between, 
or sometimes during exposures. In this mode, the num- 
ber of frames is larger, 15000 frames in our case. The cy- 
cle rate is close to 10 ms, so we recorded 15000 frames 
in about 2.5 min. We stress that MIDI is not intended 
to be an imager instrument but a long-baseline interfer- 
ometer. Therefore the majority of the targets are totally 
unresolved by a single 8m telescope, providing a wealth of 
instrumental Point Spread Function (PSF) images. The 
PSF files were recorded with the Acquisition_Chop mode 
and contain 2000 frames (20 s). 

In the following section we present the deconvolution 
treatment applied to the acquisition images. 



2.1. Images 

The data used to obtain a deconvolved image of 
01126.5-1-0.6 are summarized in Table ^ The observa- 
tions were recorded during the acquisition process and the 
source location within the field of view can be different for 
each file. The PSFs are generally well centered except for 
PSFl. Starl was very far from the FOV center, and the 
quality of the deconvolution using this observation is very 
low (but note that the results are consistent with the other 
measurements) . 

Numerous observations of two PSFs (HD 168454 and 
HD 177716) have been performed before and after the 
star acquisition. HD 168454 is a bright K3IIIa star ex- 
hibiting an IRAS 12 fiux of 62 Jy (the IRAS fiux 
of OH26. 5-1-0.6 is 360 Jy). HD 177716 is a Kllllb which 
has been observed by IRAS with a fiux of 26.9 Jy. There 
is a Cohen template available from the ISO primary cal- 
ibration database^ (Cohen et al. 1999). The visual see- 
ing during the HD 168454 exposures was ^0'.'4, during 
the OH26.5H-0.6 exposures it was ~0'.'5 and during the 
HD 177716 around 0'.'6. The airmass of the three targets 
ranges between 1 to 1.16. The pixel size on the sky is 
98mas. This scale factor has been defined from the MIDI 
observations of close visual binaries. 

The deconvolution has been performed using the Lucy- 
Richardson algorithm (1974) embedded in the IDL as- 
trolib package developed by NASA. Choosing the right 
iteration number for the Lucy-Richardson algorithm is al- 
ways a difficult task. Our goal is clearly not to perform 
the 'best' deconvolution possible but to increase the spa- 
tial resolution of the image which is well resolved by the 
UTs. The number of iterations used was between 40 and 
60. The levels where the different deconvolved images be- 
gin to disagree among each other are between 0.3% to 
1% of the maximum flux of the image, depending on the 
quality of the measurement. The level of the differences 
between PSFl to PSF6 is about 0.3%. The level of the 
differences between the PSF of HD 168454 and the ones 
from HD 177716 can reach 2% for an individual decon- 
volution but is generally 1%. PSF7 and PSFS are quite 
different, with a level of residuals reaching 1.5%. 



Table 2. Image size statistics 



Name 


FWHM X 


rms 


FWHM Y 


rms 




(mas) 


(mas) 


(mas) 


(mas) 


PSF4 


148 


24 


148 


24 


PSF6 


150 


28 


142 


20 


PSFS 


160 


26 


166 


22 


starl 


214 


8 


286 


18 


star2 


210 


2 


292 


14 


star3 


218 


4 


296 


16 


star4 


212 


2 


268 


12 



http: / / www.iso.vilspa.esa.es / users/expLlib /ISO / wwwcal / 
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Table 3. Deconvolved image parameters 



Parameter 


Mean 


RMS 


Mean radius 


240 mas 


±14 


Mean X axis 


214 mas 


±4 


Mean Y axis 


286 mas 


±6 


Mean ratio 


0.75 


±0.07 


Mean PA angle 


95° 


±6° 



The statistical properties of the PSFs have been care- 
fully inspected and a sub-set of good ones (5 over 8) have 
been selected for deconvolution. Each data set was decon- 
volved using the good PSFs and we examined the statistics 
of their geometrical characteristics. A mean deconvolved 
image was created using the mean image of star3 and star 
4 which were located at the same place in the detector. 
This reconstructed image is shown in Fig. ^ to illustrate 
the quality of the reduction process. The position of the 
target for the two other files is different and we did not at- 
tempt any shift-and-add procedure to create a single mean 
image since the useful information is extracted from a fit- 
ting procedure of the individual images. We performed a 
2D gaussian fit for each mean deconvolved image which 
provides the image position, extension and the angle of 
the long axis on the detector frame. Table 2. presents 
the statistics of this 2D gaussian procedure and Table 3. 
presents the mean parameters of the deconvolved images. 

We are confident that the star is indeed resolved at 
8.7 iim and this fact is settled definitely by looking at the 
FWHM of the spectra fSect. lT^ . However, it is necessary 
to carefully check whether the asymmetry of the image is 
real or not. The image asymmetry is strong and detected 
with a large confidence but the PA angle is almost coin- 
cident with the chopping direction. Several checks were 
performed to ensure this detection: 

1. First, as a comparison, the ratio between Y and X ex- 
tension for the PSFs is 0.997±0.05. Moreover the angle 
of the 2D Gaussian used for the fit fluctuates randomly 
and there is no indication in the PSF files that the 
chopping had any influence on the PSF's shape, i.e. 
that the chopping flagging was uncertain during the 
exposures. 

2. Second, the star is very bright and we have tested the 
deconvolution process in some carefully chosen individ- 
ual frames (4 ms exposure) taken in the middle of the 
chopping cycle. The asymmetry is already detectable 
with a SNR larger than 5 in the best quality frames. 

3. Third, we have checked in the literature whether such 
an asymmetry could have been detected in NIR by 
speckle interferometry in the past. Some asymmetries 
have indeed been reported by Mariotti et al. (1982), 
Cobb & Fix (1987) and Fix & Cobb (1988) in L, M 
and N bands. However the axis ratios detected are not 
convincing, usually within the estimated error bars of 
the measurement and not free of any bias as pointed 
out by Fix & Cobb (1988). Other speckle measure- 



0.5 



CO 

(J 

< 

-0.5 



-1 .0 

-0.5 0.0 0.5 
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Fig. 1. Contours of the mean of star 3 and star 4 decon- 
volved images. The contour levels are linearly spaced for 
the double square root of the image l^^^. The last contour 
is equivalent to 25% of the maximum of i.e 0.4% of 
the maximum of I. The three last contours are the most 
susceptible to reconstruction artifacts. The North is up 
and the east to the left. 



ments in the L and M bands are reported by Starck et 
al. (1994) based on observations carried out with the 
3.6m telescope of ESO /La Silla at pulsation phase 0.22 
(JD=2,448,429), i.e. close to maximum luminosity. A 
strong asymmetry is detected in the L band with a 
N-S/E-W ratio of the order of 0.82±0.03 after remov- 
ing the unresolved object (Starck, private communica- 
tion). It must be stressed that this measurement has 
been performed by using more than 5 orientations on 
the sky preventing any direction dependent bias. The 
agreement between their reconstructed L band image 
and our 8.7 /im image is convincing as shown in Fig. 13 
4. Fourth, surprisingly the asymmetry reported in this 
paper is correlated with the strong one reported at 
1612MHz by Bowers & Johnston (1990) at a much 
larger scale (few arcsecs). Even more interesting, 
they reported a rotation of OH26.5+0.6's shell 
at low velocity (2-3 km.s~^) for which the pro- 
jected axis is oriented along the north-south 
direction. This axis is aligned with the minor axis 
of the L and 8.7 /im images and the consequences of 
such a correlation will be discussed more extensively 
in Sect.H 

Based on the above considerations, we are convinced 
that the measured flattening is real. The data are not af- 
fected by any bias influencing the shape of the resulting 
images and this asymmetry has also been seen in other 
independent data sets. 
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Fig. 2. Contours of the reconstructed L image from the 
speckle observations at the ESO 3.6 m telescope on La 
Silla (courtesy of Starck et al.). The contours levels are 
linearly spaced for the square root of the image V-/"^. 

2.2. Spectrum 

The photometry extracted from UTl and UTS is intended 
to calibrate the recorded fringes. Two photometric files are 
recorded for each target. In the first file, one shutter only 
is opened (corresponding to UTl) and the flux is then 
split by the MIDI beam splitter and falls onto two differ- 
ent regions of the detector. The same procedure is then 
applied with UT3. The data used to get photometrically 
calibrated spectra and fringes of OII26.5+0.6 is listed in 
Table 4.. 

An independent calibration is performed for the indi- 
vidual spectra from each part of the detector and for each 
telescope. The first step is to read in the photometric data 
sets, average the frames on the target and the frames on 
the sky, and subtract the average sky frame from the av- 
erage target frame. The position of the spectrum is then 
measured column- wise by searching for peaks that are suf- 
ficiently high above the background fluctuations. The re- 
sult is the position and width of the spectrum as a function 
of wavelength. 

We use IID177716 as absolute flux calibrator (Cohen 
et al. 1999), taking into account differences in air mass be- 
tween calibrator and OH26.5. Then, the calibrated spectra 
are combined in order to provide a high SNR spectrum. 
The shape of the spectra from the same telescope agree 
within 1-2%, but the spectra from two different telescopes 
can vary by about 5%. This is due to different optical 
paths which are intrinsically different during the early use 
of the MIDI instrument with the VLTI (poor pupil trans- 
fer). This defines the limit of relative error in the shape of 
the spectrum (pixel to pixel and in terms of slope) which 
is below 5%. This limit has also been checked by extract- 
ing the spectra of several spectrophotometric calibrators 
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Fig. 3. Calibrated spectrum from MIDI (solid line) corre- 
sponding to the mean flux from UTl and UTS. The flux is 
almost 300% higher than the flux observed by ISO, shown 
with a dotted line. The ISO data have been recorded dur- 
ing a minimum of the lightcurve. 

observed by MIDI in several observing runs. The tempo- 
ral flux variations are the dominant source of error for the 
absolute flux calibration, and variations of 5-20% or even 
more are routinely observed in the N band. During the 
night of the OII26.5-f0.6 observations, the photometric 
errors were limited to 8%. 

Table 4. Journal of observations: fringes (Frg) and pho- 
tometric files (Phot) 



Star 


Telescope 


Time 


Frames 


File 


HD 168454 


UTl 


06:40:47 


3000 


Phot 


HD 168454 


UTS 


06:43:12 


3000 


Phot 


OH26.5-f0.6 


UT1/UT3 


07:13:36 


12000 


Frg 


OH26.5-f0.6 


UT1/UT3 


07:20:09 


9000 


Frg 


OH26.5-f0.6 


UT1/UT3 


07:23:46 


9000 


Frg 


OH26.54-0.6 


UTl 


07:27:51 


3000 


Phot 


OH26.5+0.6 


UT3 


07:29:51 


3000 


Phot 


HD 177716 


UTl 


08:18:44 


3000 


Phot 


HD 177716 


UTS 


08:20:48 


3000 


Phot 



We studied the spatial extension of the spectra in 
the direction of the slit in order to check if the shell of 
OII26.5+0.6 is spatially resolved at all wavelengths. A ID 
gaussian fit was performed for each column of each spec- 
trum from the target and the calibrators. The PA angle 
of the slit at 72° is close to the PA angle of the major axis 
detected in the deconvolved image at 95°. In Fig. 21 we 
see that OH26.5-I-0.6 is well resolved by the 8m telescope. 
Moreover the star is much larger in the silicate band. We 
stress here that there was no image sharpening applied, 
this result being directly extracted from the mean MIDI 
spectra. 

In order to constrain the true size of the object in the 
slit direction, we performed a deconvolution on each of the 
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Wovelength in micron 

Fig. 4. FWHM of the star spectrum from UTl (solid lines) 
throughout the wavelength range compared to the FWHM 
of the calibrator spectrum of HD 168454 (dotted lines). 
There are two lines per target because the MIDI beam 
splitter is inserted and the light falls onto two different 
regions of the detector. 




Fig. 5. Mean of FWHM curves of the star deconvolved 
spectra from UTS (solid lines) and the individual deconvo- 
lution using different calibrators (dotted lines). The errors 
bars of the figure (spaced by 0.5 /im intervals) represent 
only the scatter of the measurements. 



4 spectra available, two for each telescope. A ID deconvo- 
lution using the Lucy-Richardson algorithm is performed 
column by column using a normalized column from the 
calibrators as PSF. The same number of iterations is ap- 
plied for all the wavelengths. There are systematic differ- 
ences between the shapes provided by UTl and UTS which 
can be attributed to differences of the optical quality of 
the different light paths. 



2.3. Fringes 

Fringes were searched for by repeatedly scanning a large 
range in optical path difference between the two tele- 



scopes, using a small instrumental delay line (see Leinert 
et al. 2004 for a description of the observing sequence). 
A mask is created with the average position and width of 
the spectra recorded in the photometric files for UTl and 
UTS. This mask is used to extract the object data from 
the fringe tracking datasets. Each frame of the fringe data, 
corresponding to one individual OPD setting inside a scan, 
is reduced to a one-dimensional spectrum by multiplying 
it with the mask and performing the weighted integral 
over the direction perpendicular to the spectral disper- 
sion. Then the two - oppositely phased - interferometric 
output channels of the beam combiner are subtracted from 
each other. This combines the interferometric modulation 
of both channels into one and at the same time helps sub- 
tracting out the background. The few dozen spectra from 
each scan with the piezo-mounted mirrors are collected 
into a two-dimensional array with optical wavelength and 
OPD as axes. The contents of this array are column-wise 
fourier-transformed from OPD to fringe frequency space. 
As a rule, four of the w0.05 fim wide wavelength (pixel) 
channels were added to improve the signal-to-noise (S/N) 
ratio. The fringe amplitude for each optical wavelength 
is then obtained from the power spectrum at the corre- 
sponding fringe frequency. 

No correlated flux has been detected with the UTl- 
UTS projected basehne of 102.4 m at a PA angle of 39.6°. 
The atmospheric conditions and the data recorded dur- 
ing the fringe search have been carefully checked. The 
seeing degraded slowly between 6h and 8h UT from 0'.'4 
to 0'.'6 and the standard deviation of the flux from the 
target pointed by the (visible) seeing monitor^ increased 
also, affecting the observations of the bright calibrator 
HD177716, but still well below the cloud alert threshold. 
During all the night, the atmosphere turbulence was quite 
rapid with a mean tq—S ms. These atmospheric conditions 
while not being excellent can although be considered as 
normal conditions at Paranal observatory. Therefore the 
fringe detection threshold for MIDI during the observa- 
tions OH26. 5-1-06 was nominal. 

Based on the first few months of routine observations 
with MIDI, we can set limits on the amount of correlated 
flux the instrument is capable of detecting under average 
weather conditions. For instance, a careful data reduction 
of the data from NGC 1068 shows that a correlated flux 
can be confidently detected down to 0.5 Jy for faint objects 
(Jaffe et al. 2004). For bright objects, visibilities of the 
order of 1% have been detected from the heavily resolved 
Herbig star HD 100546 (Leinert et al. 2004), or from the 
clumpy environment of the supergiant rj Car (Chesneau 
et al. 2004). It is difficult to reach a sensitivity less than 
1% for bright objects because the beam combination is 
not perfect: a part of the noise residuals depends on the 
photometric noise from the bright source. This number 
has to be compared to the photometric flux integrated 



^ This information has been extracted from the ESO 
Ambient conditions database of Paranal observatory: 



,http://archive.eso.org/^ 
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over OH26.5+0.6 of about 600-800 Jy. With the 100m 
baseline, MIDI is sensitive to the emission coming from 
any structure smaller than 10 mas exhibiting an integrated 
flux larger than 6-8 Jy in this case. 

3. Modelling the circumstellar environment 

3.1. The approach to modeling the object 

Most of the studies on OH/IR stars and OH26.5-F0.6 in 
particular rely on the interpretation of the observed vari- 
able Spectral Energy Distribution (SED) by comparing it 
to a synthetic SED, computed using a radiative transfer 
code. In this way, one tries to separate effects of opacity 
and radial structure of the wind. Unfortunately, fits to the 
SED are often not unique, especially when deviations from 
spherical symmetry are taken into accout. MIDI however 
provides a unique spectrally and spatially resolved data 
set which puts strong new constaints on any model for 
the envelope of OII26.5+0.6. Below, we briefly outline the 
strategy used to fit the MIDI observations. 

Despite the convincing evidence that the envelope of 
OH26.5-I-0.6 is not spherical, we begin our analysis assum- 
ing spherical symmetry. We use the SED observed by ISO 
to determine global envelope parameters. The ISO spec- 
trum (and unfortunately also the IRAS data) was taken 
close to the minimum luminosity of the star (Suh & Kim 
2002). Therefore the model fit provides some constraints 
on the physical parameters of the dust shell close to the 
minimum luminosity. The model fit parameters are com- 
pared to the parameters published by Justannont et al. 
(1996), who also mainly scaled their spectrophotometric 
data to the minimum phase. 

As a second step, we tried to find a good fit to the MIDI 
data alone (taken at maximum luminosity), i.e. without 
the help of any external spectrophotometric information 
as performed by Suh & Kim (2002). This is, we tried to fit 
the MIDI spectrum by performing slight modifications to 
the minimum light model. Our goal was to check whether 
the MIDI spectrum can be fitted based on the previous 
model. As soon as a satisfactory fit is reached for the MIDI 
spectrum we evaluate the spatial distribution of the flux 
predicted by the model and compare it with the extension 
at each wavelength measured by MIDI. 

Clearly it is not is the scope of the paper to find the 
best model for the dust shell at maximum luminosity and 
such a complex study is let to a dedicated paper. Our 
goal is to pinpoint the kind of constraints provided by the 
inclusion of the spectrum and the spatial extension of the 
object in the process of model fitting and demonstrate how 
new information can emerge on the dust content close to 
the star. 

We use the radiative transfer code, MODUST, com- 
monly used for the SED fitting of this kind of stars. The 
radiative transfer technique implemented in this code has 
been outlined by Bouwman (2001) and the specification 
of grain properties, such as size and shape distribution, is 
discussed in Bouwman et al. (2000). 
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Fig. 6. Comparison between the ISO spectra (SWS + 
LWS) and a MODUST model using the JU96 parameters 
but with a slightly reduced luminosity of the central star. 
Other differences arc the use of CDE theory and the inclu- 
sion of metallic iron. The crystalline features present in the 
ISO-SWS observation are not included in the model since 
they have little effect on the model structure (Kemper et 
al. 2002) 



Throughout this analysis, we will adopt a distance of 
1.37 kpc for OH26.5-I-0.6. We note that distances to AGB 
stars are notoriously uncertain, and OH26.5-I-0.6 is no ex- 
ception. 

3.2. SED fit of OH26.5+0.6 at minimum luminosity 

JU96 constructed a model for OH26.5-I-0.6 based on a 
large collection of photometry, spectroscopy and CO-line 
measurements. Their 2-component model (a very thick in- 
ner shell, due to a recent superwind, surrounded by a ten- 
uous AGB wind) shows that the near-mid IR is dominated 
by the superwind region. Most of their study is based on 
fluxes measured at minimum light, and hence should also 
be applicable to the ISO-SWS spectrum, also obtained 
around minimum light (JD=2,450,368, phase 0.47) 

The superwind hypothesis is confirmed by the ISO- 
SWS spectrum, which shows very little far IR flux w.r.t. 
the depth of the 10 fxm feature. Assuming a density distri- 
bution going as r~^ this can be modelled only by cutting 
the shell fairly close to the star (at a few hundred stellar 
radii instead of a few thousand). 

Looking at the shape of the 9.7 jim feature, we can al- 
ready make some improvements on the composition of the 
dust. From its width and the location of the minimum, we 
conclude that CDE theory (CDE, Continuous Distribution 
of Ellipsoids) is to be preferred over spherical dust parti- 
cles. Furthermore, there is also strong evidence for the 
presence of metallic iron, as is the case for OH127. 8-1-0.0 
(Kemper et al. 2002): the slope in the (ncar-)IR (4-8 /xm) 
cannot be explained without it. It is worth noting that the 
flux blocked by the metallic iron in the near-IR emerges 
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Table 5. Model parameters from Justtanont et al. 96 
(JU96), and the same model adapted to minimum light 
(the ISO data, JD=2,450,368, phase 0.47) and maximum 
light (the MIDI data, JD=2,452,804, phase 0.06), dust 
evaporation scenario). 



Par am. 


JU96 


Minimum 


Maximum 


Te// (K) 


2200 


2200 


2100 


(R©) 


862 


650 


1100 


Dist. (kpc) 


1.37 


1.37 


1.37 


Super wind 


Rin (R*) 


7.5 


7.5 


20 


Rout (cm) 


8 X 10^^ 


8 X 10^^ 


8 X 10^^ 


M (Me/yr) 


5.5 X 10"'' 


5.5 X 10'* 


8.5 X lO""* 




AGB 


wind 




Rin (cm) 


8 X 10'^ 


8 X 10^^ 


8 X lO^'' 


Rout (cm) 


5 X 10^® 


5 X 10^** 


5 X 10^** 


M (Me/yr) 


1. X 10-'^ 


1.4 X 10'^ 


1.4 X 10~^ 



again in the Mid to Far-IR. Hence, the amount of Fc will 
significantly influence the optimum value of the other shell 
parameters. 

The luminosity used by JU96 is by far too high for the 
epoch of the ISO-SWS observations. However, we do ob- 
tain a satisfying fit by reducing the stellar radius to 650 
Rq and keeping the outer radius at 8x10^^ cm. The com- 
parison model vs. ISO-SWS spectrum is shown in Fig.lHl 
The general shape is approximately good. Most of the dis- 
crepancies can be attributed to the lack of crystalline dust 
in our model. The crystalline features present in the ISO- 
SWS observation are probably due to a few percent of 
enstatite and forsterite but since these do not influence 
significantly the model structure (their opacities are very 
similar to those of the amorphous material, Kemper et al. 
2002), we do no detailed fitting of their spectral features. 

3.3. Attempts to account for the MIDI data 

Our model will have to explain the following new MIDI 
observations: 

1. The MIDI N band spectrum taken close to the maxi- 
mum hght of OH26.5-f 0.6 with the 0'.'6x2" slit. 

2. The spatial extent of the spectrum. We limit ourselves 
to the comparison of the FWHM provided by a fit of 
the PSF-deconvolved MIDI spectrum by a ID gaus- 
sian compared with a similar fit of the model intensity 
profiles. The intensity distribution of OH26.5+0.6 on 
the sky is likely more complex than a simple gaussian. 
However, it turned out to be difficult to disentangle im- 
perfections in the imaging quality resulting from the 
many reflections in the VLTI optical train from those 
intrinsic to the source, especially at lower intensity lev- 
els. The observed FWHM is attributed to a spherical 
object. The slit was oriented at PA=72°, i.e. close to 
the maximum extension of the object. 
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Fig. 7. Comparison between the MIDI spectrum (origi- 
nal: solid line, dereddened for IS extinction: dashed line) 
and the spectrum resulting from our ISO-tuned model but 
with increased central star luminosity (diamonds) 
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Fig. 8. A comparison between the FWHMs of the inten- 
sity profiles coming from the ISO-tuned model with a 
more luminous central star and the MIDI FWHMs. The 
predicted variation of diameter with wavelength is much 
larger than the one observed. Mainly the size at minimum 
optical thickness of the shell (at 8.5 and 13.5 jim) does 
not agree. 



3. The negative detection of fringes by MIDI. The con- 
straints provided by this observational fact should not 
be underestimated. During the maximum luminosity, 
the dust shell opacity is at its minimum. It is quite 
difficult from a model of the MIDI spectrum alone to 
disentangle models which are almost optically thin in 
the wings of the silicate features to the ones more op- 
tically thick. The MIDI observations definitely discard 
any models of OH26.5-I-0.6 for which the central star 
is visible with a (correlated) flux larger than 3-6 Jy 
through the shell at any wavelength located between 
8 jjLm. and 13.5 /xm. 
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The MIDI observations were done very close to max- 
imum light (JD=2,452,804, phase 0.06), resulting in an 
observed flux which is more than twice as high compared 
to minimum light. In order to fit the model at minimum 
luminosity to the MIDI data at maximum luminosity, we 
must increase the total luminosity to 1.7 x W^Lq. 

The increase in total luminosity is simulated by an 
increase of the central star diameter, still keeping the ab- 
solute value of the outer radius of the superwind fixed. 
Below we confront this model with the MIDI FWHM ob- 
servations. 

Fig- El shows the spatial intensity profiles according to 
our ISO-tuned model adapted to the higher total luminos- 
ity at the time of the MIDI observations. The profiles at 
8.5 and 13.5 /xm correspond to a fairly low optical thick- 
ness of the shell and thus the central star is not totally 
obscured. However, the amount of correlated flux by the 
central star is at most a few Jansky and hence close to 
the detection hmit of MIDI (1 % or 5-8 Jy). At 10.5 /im, 
the shell reaches an optical thickness of more than 10, 
resulting in the gaussian intensity profile. 

The FWHM's determined from these intensity profiles 
range from 100 mas to 370 mas (Fig. ISJ, clearly show- 
ing that if only the opacity by amorphous olivines were 
to determine the diameter seen, the variations with wave- 
length would be much larger than what is observed. The 
maximum size appears to compare reasonably well with 
the model (though with a slightly different wavelength of 
maximum) and hence is compatible with the superwind 
size of JU96. More precisely, we can put a lower limit to 
the size of the superwind region of 400 mas, which corre- 
sponds to 4 X 10^^ cm at 1.37kpc. 

The discrepancy between the spatial extent of the 
model and the observed size of the dust shell near 8 and 
13 /im can be resolved if we move the inner radius of the 
dust shell to larger distance from the star. The minimum 
size seen will be mainly determined by the inner radius 
of the dust shell: at the wavelengths of low optical depth, 
the intensity profile is not at all gaussian fFig. 110(1 . In this 
way, the observed diameter in the wings of the 9.7 /im 
feature can be simulated by moving the inner radius out 
to about 20 R0 fFig.lll|l. For this model, a slightly higher 
mass loss (8.5 x 10""* M0) is needed to keep the quality of 
the spectral fit. 

At first sight, an alternative solution would be to in- 
troduce a source of opacity with only a modest wavelength 
dependence, which would dominate over the silicate dust 
opacity near 8 and 13 /im. This could either be gas-phase 
molecular opacity or dust. However, by doing so the spec- 
tral fits become unacceptably poor, because the depth of 
the silicate feature can no longer be reproduced. 

In summary, we conclude that the spectral and spatial 
data of OH26.5-f 0.6 can be understood in the framework 
of a spherically symmetric shell, with dust components 
that are also shown to be present in other OH/IR stars. 
The outer radius of the dust shell agrees well with the 
one estimated by Justannont ct al. (1996) as the radius of 
the onset of the superwind. The spatial data near 8 and 13 
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Fig. 9. Normalized intensity profiles for our model at max- 
imum luminosity. In the wings of the 9.7 /im profile (at 
8.5 and 13.5 /im ), the shell optical thickness is only about 
2 and thus the central star is visible. At 10.5 /im, the shell 
reaches an optical thickness of more than 10, resulting in 
the gaussian intensity profile. 
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Fig. 10. Normalized intensity profiles for our model at 
maximum luminosity with an increased inner radius. For 
such intensity profiles, the inner radius is determining the 
observed size of the object at wavelengths of low opacity. 

/im force us, in the context of spherical symmetry, to move 
the inner radius of the dust shell to a distance of about 20 
stellar radii. However, this results in a predicted correlated 
flux using the 102 meter UT1-UT3 baseline which is 5 
times the upper limit imposed by the non-detection of 
fringes in the intcrfcromctric signal (under the assumption 
that the central star has a typical AGB temperature of 
3000-4000 K). 

4. Discussion 

4.1. Nature of the large inner dust radius 

While a large inner radius of the dust shell seems a sim- 
ple solution to our fitting problems, it is clearly not in 
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Fig. 11. The large observed radius at the red and blue 
sides of the profile can be simulated by moving the inner 
radius of the model far out, to about 20 R*. 

agreement with the limits set by the interferometric mea- 
surement. In addition, such a large inner radius is not 
compatible with our current understanding of oxygen-rich 
AGB dust shells: the dust temperature at the inner edge 
of our dust envelope is only 500-600 K, well below what 
is believed to be the condensation temperature of olivines 
(1000 K). Furthermore, one can wonder whether the re- 
gion between photosphere and olivine dust shell contains 
other material (refractory residuals of dust for instance, 
like corundum). Given the constraint that whatever fills 
this region must be quite transparent from 8 to 13 /xm^, 
several hypotheses can be formulated: 

— A large cavity within this inner region, could indi- 
cate that the mass loss has decreased strongly about 
30 years ago. This is compatible with the timescales 
derived from the rings observed around several Post- 
AGB stars, hinting at episodic mass-loss with peri- 
ods of a few hundred years (e.g.IRAS 17150-3224: 
Kwok, 1998;IRC+10216: Mauron & Huggins 1999; Egg 
Nebula:Sahai, 1998, Marengo, 2001). However, stellar 
pulsations are quite regular and have been detected 
over the past 30 years (e.g. Suh and Kim, 2002). So, if 
mass loss stopped, apparently the pulsations did not. 
Moreover, this model offers little opacity at 8 micron 
{TSfim is of the order of 1), and hence the amount of 
correlated flux from the central object would be of the 
order of 5 times the detection limit of MIDI. This puts 
into question the relation between pulsations and mass 
loss for AGB stars in the superwind phase, 

— The mass loss is continuing even today and no con- 
densation of the dust is possible inwards of 20 stellar 
radii when the star is at maximum luminosity. In fact, 
early dust condensation models by Sedlmayr (1989) 

^ This means optically thin from 7.5 to 9 micron and from 
11.5 to 13.5 micron. An optically thick shell (t ~2-3) would 
ruin the diameter profile since a smaller radius would be nec- 
essary to account for the SED. 



have predicted that the dust condenses only when the 
gas is extremely super-saturated, which happens well 
below the glass temperature, and would be around 800- 
600K. Several other studies have previously hinted at 
the possibility that dust formation in AGBs does not 
happen close to the star. Danchi ct al. (1994) also 
find some examples of stars with rather detached dust 
shells, corresponding to timescales of decades, so sim- 
ilar to what we find. However, other stars have inner 
dust radii much closer to the star. 
— The dust gets periodically destroyed (dust evapora- 
tion) because the stellar luminosity changes during 
a pulsation cycle. Some calculations were done by 
Suh and Kim (2002), but do not predict the required 
amount of dust destruction to agree with our geometri- 
cal model. It is not certain that the net amount of dust 
created through an entire cycle is stable. This means 
that the episodical mass loss is a matter of balancing 
between destruction and formation. Of course, if the 
dust formation is not large enough to compensate for 
the evaporation, the inner gap would increase in size 
over multiple cycles. 

It is not straightforward to choose between these pos- 
sible explanations, and more importantly, not one of them 
explains the non-detection of fringes. 

Although the SED is compatible with our spherically 
symmetric model, this assumption might be strongly vi- 
olated, as suggested by the acquisition image. It is intu- 
itively clear that a disk-|-bipolar outflow structure could 
account for the fitting problems we experienced. If the 
MIDI slit was oriented perpendicular to a nearly edge- 
on disk, the very high density (and thus opacity) would 
explain both the large size in the wings of the silicate 
feature and the non-detection of fringes. Because of the 
much larger complexity of 2D CSE modelling, we leave an 
in-depth analysis of this hypothesis to a future dedicated 
paper. 

4.2. Star asymmetry 

The drastic changes obscrvcxl in AGB and post-AGB 
stars circumstellar structure are particularly puzzling. 
Information on the geometry of the AGB mass loss was 
first obtained from interferometry maps of OH 1612 Hz 
maser emission (Booth et al. 1981; Herman et al. 1985, 
Bowers & Johnston 1990). These, as well as CO radio line 
maps of a large number of AGB-stars (Neri et al. 1998, 
Olofsson et al. 1999 and references therein), were con- 
sistent with an over-all spherically symmetric mass loss. 
Evidences concerning deviation of the geometry of OH/IR 
stars can be found in the literature but the results are 
generally somehow contradictory and perplexing due to 
the large range of masses and evolutionary stages encom- 
passed by this term, from embedded AGBs to PPNs (Van 
Winkel, 2003). In particular, the sub-group of OH/IR star 
for which no pulsation can be detected is associated with 
PPNs (Zijlstra et al. 1991). 
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What makes the MIDI observations particularly inter- 
esting is that JU96 demonstrated that the star entered 
the supcrwind phase ^200yr ago. This provides an ex- 
tremely short upper limit to the development of a large 
scale asymmetry. For a large mass star like 0H26.5-I-0.6, 
the duration of the OH/IR phase is expected to be of the 
order of lO^'yr, to be compared to lO'^yr for a low-mass star 
(JU96). Increasing evidence show that AGB wind becomes 
axi-symmctric at the very last stages of the AGB evolu- 
tion, and that the interaction with a fast wind from the 
post-AGB object further enhances the axi-symmetry. The 
MIDI observations suggest that in the case of OH26.5-I-0.6 
the appearance of asymmetries can occur is a fairly short 
time scale (Sahai et al. 2003). 

JU96 and Fong, Justtanont, Meixner, Campbell (2002) 
reported CO observations which did not show any signifi- 
cant deviation from spherical symmetry for OH26.5-I-0.6, 
but at most of the emission is spatially unresolved (com- 
ing mostly from the superwind). In contrary, the OH 
1612 Hz maser emission from Baud (1981) and Bowers & 
Johnston (1990) present a clear picture of the clumpy and 
asymmetric environment OH26.5-I-0.6. The radio shell of 
OH26.5+0.6 is certainly one the most extended and least 
symmetric known for OH/IR stars. The large scale envi- 
ronment of OH26.5-I-0.6 is crowded which not ease the ex- 
traction of the radio emission sensitive to any anisotropic 
UV radiation field and Bowers & Johnston (1990) pro- 
posed that it is a likely cause for the detected asymmetry. 

The crucial point is that the axis of symmetry of the 
present Mid-IR objet is perfectly aligned to the large scale 
anisotropy detected in radio, which excludes a-priori a 
strong external influence on the shaping of the OH maser. 
This correlation opens new possibilities of interpretation. 
Hence, Bowers & Johnston (1990) detected some hints of 
rotation at low projected velocity (v^ <3 km s— 1) with a 
rotational axis aligned with the minor axis of the asym- 
metric shell. 

What could be the origin of the development of such 
an asymmetry? It is not in the scope of this paper to 
review all the mechanisms invoked for explaining such a 
phenomenon and the reader is invited to consult the re- 
view from Balick & Frank (2002). We just point out that 
OH26.5-I-0.6 is not a known binary, but owing to the dif- 
ficulties to study the central star of OH/IR star, this lack 
of detection is not meaningful. The (iisc;ussion above lead 
us to think that OH26.5-I-0.6 particular characteristics are 
perhaps better understood under the binarity hypothesis. 

4.3. Improving the model 

By using an up-to-date spherical model of a dust shell, we 
have been able to fit satisfactorily the SED of the star but 
this model failed to provide a direct explanation of the 
non-detection of any fringes within the N band. 

These difficulties point to a problem of opacities lo- 
cated in regions fairly close to the star though sufficiently 
extended to prevent the detection of correlated flux by 



MIDI. One of the remedy is to make the radius of the 
central source larger so that almost no correlated fluxes 
can be detected by MIDI with a baseline as large as 100m. 

An attractive solution to this problem would be the in- 
clusion of molecular opacities. Growing evidence of their 
deep effects on interferometric measurements in the near 
and mid-infrared are reported (Matsuura et al. 2002, 
Mennesson et al. 2002, Perrin et al. 2004a, b, SchuUer 
et al. 2004, Cotton et al. 2004, Ohnaka et al. 2004). The 
first effect is to increase the diameter of the central star 
and thus, to decrease the correlated flux. At maximum 
luminosity, the expected angular diameter OH26.5-I-0.6 is 
about 8 mas (for R*=11OOR0) and the correlated flux 
from the central object should represent about 80% of the 
stellar flux if the star is a uniform disk. The inclusion of 
an optically thick molecular envelope of H2O and SiO of 
about 2.5R* divides this correlated flux by 10, prevent- 
ing probably its detection by MIDI. Moreover, the star 
can probably no longer be modelled by a uniform disk 
but by a spatially smoother flux distribution which de- 
creases again the correlated flux. The second effect is to 
redistribute the flux from the central star to other regions 
of the spectrum. Of course, if the warm molecular layers 
are optically thick they will emit like a blackbody at a 
temperature slightly lower than the star. The effects of 
this envelope on the dust formation/destruction processes 
have to be carefully evaluated and need consistent radia- 
tive transfer calculations which are not in the scope of this 
paper. 

The opacity of the dust shell remains an issue which 
could be solved by changing the dust density without af- 
fecting its composition or by adding species particularly 
absorptive between 8 and 13.5/im. This implies either a 
very large olivine dust column density, cither other highly 
absorptive dust species like corundum. An inner shell rich 
in metallic iron and corundum (AI2O3) could help to fit 
the observed extensions, preventing the flux from the star 
to be detected. The corundum is used for instance to 
model the opacities of the thin dust shell around some 
Miras, with mass ratio of cormidmii to silicate ranging 
from 0.6 to 0.9 and grain sizes of 0.1 - 0.2 fim (Martin- 
Lorenz & Pompeia, 2000, Ohnaka et al. 2004). However, 
the amount of material needed is limited by the natural 
abundance of Al in the photosphere. In our first tests, 
this amount needed is at the moment unrealistically high 
to prevent the correlated flux from the star to be detected, 
under the hypothesis that the star is a naked photosphere 
of about IIOOR0 at 2100K (i.e. a uniform disk). 

Finally, for the sake of simplicity, we have put all the 
discussions in the frame of a spherical object. All the 
codes used to model the SED of OH26.5-I-0.6 have been 
using spherical geometry to understand an object which 
is proven to be strongly fiattened in this article. A very 
promising hypothesis is that we are indeed observing a high 
olivine column density in the direction perpendicular on 
the slit, i.e. that we are looking an equatorial overdensity 
(or even a disk) close to an edge-on configuration, explain- 
ing the important aspect ratio of the 8. 7fjim MIDI image. 



12 



O. Chesneau et al.: The Mid-IR spatially resolved environment of OH26.5+0.6 at maximum luminosity. 



4.4. Envelope clumping 

The fact tliat no fringes have been detected from 
OH26.5+0.6 implies also that the dusty environment of 
OH26.5+0.6 is relatively homogeneous and smooth. Most 
of the flux originates from the dust shell and the absence 
of fringes is in great contrast with the almost ubiquitous 
fringes found around the massive star Eta Car in an area 
as large as 0'.'6x0'.'6 with photometric fluxes comparable 
to the ones reported in this paper (Chesneau et al. 2004). 
The use of continuous dust distribution for the modelling 
of this kind of environment is thus fully justified. 

The pulsations are supposed to generate a strongly 
clumped medium due to the shocks, but Suh et al. (1990) 
has shown that this region is limited to the 3 R* . The rapid 
outward acceleration extending to 10-20R* should consid- 
erable smooth the dusty wind. In the Suh et al. model, 
the dust condensation radius is about 6 R* depending on 
the pulsation phase. Even considering the pulsation, their 
model of the dusty envelope is very close to a smooth r^^ 
density law. The clumpy regions embedded in the optically 
thick part of the shell at 10 /um should not be visible but 
some signal could be expected at 8 or 12 /im at maximum 
luminosity if the dust shell is sufficiently optically thin. We 
note that the clumps are probably embedded in the puta- 
tive optically thick molecular layer so that their correlated 
emission would largely hidden and therefore undetectable 
by MIDI. 

5. Conclusion 

It has been shown that the dust model used in this article 
to interpret the SED from 01126.5+0.5 has strong difficul- 
ties to predict the large extension of the dust shell outside 
the silicate absorption region with simultaneously having 
opacities sufficient to render the flux from the central ob- 
ject undetectable by MIDI. The conjunction of these two 
complementary constraints make it necessary to pursue a 
deeper effort to understand the physical processes operat- 
ing at the inner regions of the dust shell. 

OH26.5+0.6 is indeed a very complex object exhibiting 
a wide range of physical phenomenon: 

— an asymmetric appearance, whose axis of symmetry 
is probably coincident with the axis of rotation of the 

star, 

— a thick dust envelope whose characteristics are modu- 
lated by the pulsation cycle, 

— a complex inner shell where dust forms and is destruc- 
ted throughout the cycle. The contribution of corun- 
dum and metallic iron opacities in this region are prob- 
ably important, 

— a putative thick molecular envelope as encountered in 
many, if not all, Mira stars which increases the an- 
gular diameter of the central star, and decreases the 
apparent temperature, 

The results presented in this paper are very constrain- 
ing and have to be confronted with a model able to handle 



consistently the complex interplay between the pulsating 
central star, its molecular atmosphere and the mechanisms 
of dust formation/destruction and transport. However, 
such a theoretical approach is at the moment inefficient 
until any confirmation of the presence of a disk around 
this object allows a restricted range of geometrical pa- 
rameters. 

In the course of the ^1560 days of the cycle, 
MIDI/VLTI interferometer will allow a continuing mon- 
itoring of OH26.5-I-0.6. The observations will provide a 
unique view of the evolution of the size and the shape 
of the dusty envelope throughout the entire cycle. However 
the phase of maximum luminosity remains the unique op- 
portunity to reach the most internal regions close to the 
star. 

High resolution observations in optical and near- 
infrared by means of Adaptive Optics should also help to 

estimate the amount of scattered light close to the object 
in order to test efficiently the disk hypothesis. 
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